Phase transitions and iron-ordered moment form factor in LaFeAsO 
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Elastic neutron scattering studies of an optimized LaFeAsO single crystal reveal that upon cooling, 
an onset of the tetragonal (T)-to-orthorhombic (O) structural transition occurs at Ts w 156 K, and 
it exhibits a sharp transition at T P « 148 K. We argue that in the temperature range Ts to Tp, T and 
O structures may dynamically coexist possibly due to nematic spin correlations recently proposed 
for the iron pnictides, and we attribute T P to the formation of long-range O domains from the finite 
local precursors. The antiferromagnetic structure emerges at T N ~ 140 K, with the iron moment 
direction along the O a axis. We extract the iron magnetic form factor and use the tabulated (jo) 
of Fe, Fe 2+ and Fe 3+ to obtain a magnetic moment size of ~0.8 /xb at 9.5 K. 

PACS numbers: 74.25.Ha, 74.70.Xa, 75.30.Fv, 75.50.Ee 



I. INTRODUCTION 

The discovery of _RFeAs(Oi_ x F x ) (R = rare earth, 
"1111") superconductors pQ with transition tempera- 
tures up to 56 K [21 [3] has stimulated a renewed ex- 
citement in the search for novel superconductors derived 
from antiferromagnetic (AFM) parent compounds. In 
these iron-arsenide based compounds, the superconduct- 
ing (SC) state can be achieved by doping or by the appli- 
cation of pressure. The appearance of superconductivity 
is normally accompanied by a suppression of the AFM 
state in the parent compounds and by the appearance of 
a spin resonance [H[S]. It is thus important to understand 
the magnetism in these iron-based parent compounds in 
order to unravel the mechanism that leads to supercon- 
ductivity. 

Because of the initial difficulties to grow large "1111" 
single crystals, most of the research on these systems has 
been performed on polycrystalline samples. These stud- 
ies showed that LaFeAsO undergoes a structural phase 
transition from the tetragonal (T, PA/nmm) symmetry 
[Fig. flja)] to the orthorhombic (O, Cmraa) one [Fig. 
[lib)] at T s « 155 K and forms an AFM ordering (also 
referred to as a spin-density-wave - SDW) at T N f=a 137 
K upon cooling, with an ordered magnetic moment of 
0.36(5) /!b [6J. Implied in these reports on powdered poly- 
crystalline samples [6-8J is that the form factor used is 
that of Fe 2+ ion. However, first-principle calculations on 
LaFeAsO predict a larger localized magnetic moment of 
^2.6 /is at each iron site that is embedded in an itiner- 
ant electronic environment [§]. It should be noted that 
density functional theory argues that a large enough mag- 
netic moment ~2 /u B is necessary to drive the observed 
O-T transition [TUj. It is clear that there is a significant 
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(a) PA/nmm 



FIG. 1: (color online) (a) Crystal structure (P4/nmm) above 
T s = 156(1) K with two unit cells (solid lines), (b) Crystal 
structure (Cmma) below T s with two unit cells (solid lines) 
and AFM structure below T« = 140(1) K in one AFM unit 
cell. The arrows on the Fe irons in (b) represent the spins 
of iron in single-crystal LaFeAsO. The unit cells of P4/nmm 
(T), Cmma (O) and AFM structures are (a b c) (a — b), 
(y/2a y/2b c) and (y/2a y/2b 2c), respectively. 



discrepancy between the calculated iron moment and the 
experimental one. Therefore, determining the moment 
size is vital to the validity of any theoretical model that 
attempts to explain the electronic structure of LaFeAsO. 
Another important issue is the coupling between struc- 
tural and magnetic behaviors in iron pnictides. Inelastic 
neutron scattering studies from polycrstalline LaFeAsO 
showed two-dimensional magnetic fluctuations that per- 
sist up to room temperature (over ^160 K above T N ). It 
was argued [TlTfTS"] that such fluctuations introduce dy- 
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FIG. 2: (color online) Observed (dots) and calculated 
(solid line) X-ray powder diffraction patterns for pulverized 
LaFeAsO single crystals at ambient conditions obtained on an 
in-house diffractometer employing the copper K a i = 1.54056 
A and K a 2 — 1.54439 A with I a 2/Iai — 0.5 as the radiation. 
The vertical bars mark the positions of Bragg reflections. The 
lower curve represents the difference between the observed 
and calculated patterns. 



namic disorder of the O/T mixed-phase in the so-called 
T phase, suggesting that a finite orthorhombicity may 
exist above the O-T structural transition [TlTU3| . 



TABLE I: Refined structural parameters with T (P4/nmm) 

symmetry for pulverized LaFeAsO single crystals at room 

temperature. 

a = 4.0316(1) A, c = 8.7541(1) A, V = 142.290(1) A 3 

R B = 1.01, R F = 0.60, x 2 = 1-18 



atom 


site 
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B(k 2 ) 


La 


2c 


0.25 


0.25 


0.1405(2) 


1.65(4) 


Fe 


2b 


0.75 


0.25 


0.5 


1.69(12) 


As 


2c 


0.25 


0.25 


0.6543(3) 


2.15(8) 


O 


2a 


0.75 


0.25 


0.0 


insensitive 



The polycrystalline studies have provided important 
insights into the behavior of LaFeAsO, but questions re- 
main with regard to the AFM moment direction, the 
moment size and its spatial distribution. To a large ex- 
tent these questions are related to the complexity of iron 
chemistry, namely its valence (or mixed- valence) , bond- 
ing, and electron configuration in different local chemical 
surroundings. To address these questions, the study of 
high-quality single crystals is vital. Such crystals are now 
available with the recent successful growth of relatively 
large "1111" single crystals at ambient pressure [2]. 

Here we report elastic neutron scattering and syn- 
chrotron X-ray powder diffraction studies on the bulk and 
pulverized LaFeAsO single crystals, respectively, focusing 
on the details of the structural and magnetic transitions 



as well as the coupling between them, and the measure- 
ment of the average ordered magnetic moment size and 
its spatial distribution by extracting the magnetic form 
factor of iron in this compound. 



II. EXPERIMENTAL 

LaFeAsO single crystals were synthesized in an NaAs 
flux at ambient pressure as described in a recent re- 
port [2]. Crystal quality was characterized by Laue 
back-scattering, X-ray powder diffraction, heat capac- 
ity, magnetization and resistivity measurements. A large 
LaFeAsO single crystal (~20 mg) was selected for this 
study. The mosaic of this single crystal is 0.59(2)° full 
width at half maximum (FWHM) for the (202) (O nota- 
tion) reflection at room temperature. The elastic neutron 
scattering measurements were carried out on the HB-1A 
fixed-incident-energy (14.6 meV) triple-axis spectrome- 
ter using a double pyrolytic graphite (PG) monochroma- 
tor (located at the High Flux Isotope Reactor, HFIR, at 
the Oak Ridge National Laboratory, USA). Two highly 
oriented PG filters, one after each monochromator, were 
used to reduce the A/2 contamination. The beam colli- 
mation throughout the experiment was kept at 48'-48'- 
sample-60'-360'. The single crystal was wrapped in Al 
foil and sealed in a He-filled Al can which was then loaded 
on the cold tip of a closed cycle refrigerator with (h0l) 
in the Cmma symmetry as the scattering plane. The 
synchrotron X-ray powder diffraction study of pulverized 
LaFeAsO single crystals from the same batch as the one 
used for the neutron diffraction study was carried out at 
the 11-BM beamline, Advanced Photon Source, Argonne 
National Laboratory. 11-BM is a bending magnet beam- 
line, equipped with a vertical beam collimation mirror, 
a double crystal monochromator with a horizontal sagit- 
tal focusing second crystal, and a vertical focusing mir- 
ror. The calibrated X-ray wavelength was 0.41219(1) A. 
We note that the (HKL)-? indices for T symmetry corre- 
spond to the O reflection (hkl) based on the relations 
h = H + K, k = H - K, and I = L. 



III. RESULTS AND DISCUSSION 
A. Structural and magnetic transitions 

A room-temperature X-ray powder diffraction pattern 
of pulverized LaFeAsO single crystals from the same 
batch as the one used in this study measured on an in- 
house diffractometer is shown in Fig. [2] (dots), where the 
best structure refinement with P4/nmm symmetry using 
Fullprof suite (15] is also displayed (solid line) , yielding a 
good fit (x 2 = 1.18) and indicating a high degree of phase 
purity of the single crystals. The refined structural pa- 
rameters listed in Tab. Ill are in agreement with previous 
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FIG. 3: (color online) Temperature dependence of the nor- 
malized integrated intensities of (a) 8-28 (void symbols) and 
rocking-curve 8 (solid symbols) neutron diffraction scans of 
the nuclear Bragg (400)o (circles) and (006)o (triangles) re- 
flections, and (b) the rocking-curve scan of the AFM (103) 
neutron diffraction reflection (down-solid-triangles) of single- 
crystal LaFeAsO measured upon warming up the crystal from 
9.5 K. (c) and (d) are the enlargements of (a) and (b) near the 
structural and AFM transitions, respectively. The integrated 
intensities of (400)o and (006)o rocking curves are rescaled by 
1.66 and 2.10, respectively, as indicated to make them coin- 
cide with the respective 8-28 scans at low temperatures. 



reports [71 Hi]. 

Figure [3] shows the temperature evolution of the struc- 
tural and magnetic transitions in the single crystal mea- 
sured by neutron diffraction on HB-lA. For the struc- 
tural transition the (400) /(040) and (006) D (Cmma) 
reflections were monitored. The integrated intensity of 
the rocking-curve scans of the (400) reflection shows an 
appreciable increase at — 156 K, with a relatively sharp 
maximum at —148 K. The scattering at the (400) /(040) 
rocking curves or 9-28 scans is strongly influenced by 
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FIG. 4: (color online) Synchrotron X-ray powder diffraction 
data of pulverized LaFeAsO single crystals above (158 K), 
during (150 K) and below (130 K) the O-T transition (circles). 
The solid lines are fits of the Lorentzian line-shape. 
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FIG. 5: (color online) Temperature variation of FWHM 
of the (004)o/t (a) and (220)t (circles) (b) reflections from 
synchrotron X-ray powder diffraction study of pulverized 
LaFeAsO single crystals. Below the O-T structural transi- 
tion, the FWHM of the (220) T reflection was averaged from 
the (400)o and (040)o reflections. To avoid the influence of 
systematic errors, we normalized the FWHM of the (220)t 
reflection to that of the (004) /t as shown the squares in (b). 
The lines in (b) are guides to the eyes. 



TABLE II: The q values, integrated intensities, and the cor- 
responding values of \F K (hkl)\ 2 sin 2 a at 9.5 K of the AFM 
reflections observed in single-crystal LaFeAsO with the AFM 
structure as shown in Fig. fTfb). 

\F H (hkl)\ 2 sin 2 a 
6.138 
31.251 
46.470 
53.670 
57.325 
6.138 
14.560 
59.373 
23.424 



TABLE III: The fitted results of Fig. [6] with different form 
factors (jo) of Fe, Fe 2+ and Fe 3+ with spin contribution only 
[16] . and a homogeneous spherical model (details in text): the 
corresponding magnetic moment size M (/jb) and ionic radius 
R (A). The effective ionic radius of Fe 2+ and Fe 3+ was taken 
from 22^ with the CN = 4. 
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q (A" 1 ) 


Intensity 


(101) 


1.162 


26.55±0.35 


(103) 


1.545 


107.50±1.00 


(105) 


2.111 


94.08±0.81 


(107) 


2.751 


53.87±0.48 


(109) 


3.422 


22.21±0.31 


(303) 


3.487 


2.98±0.11 


(305) 


3.773 


3.35±0.11 


(1011) 


4.110 


7.29±0.13 


(307) 


4.164 


2.42±0.08 



Models 



M (us 



R(A) 



SrFe2As2/bcc-iron 

Fe 2 + 

Fe 3 + 

sphere 



0.83±0.04 
0.82±0.03 
0.77±0.03 
0.74±0.04 



0.63 

0.49 

0.63±0.04 



twinning formation upon cooling [13] that may increase 
the in-plane mosaicity in the crystal due to nucleation 
of domains with different orientations at different sites. 
Based on these observations, we suggest two critical tem- 
peratures for the structural transition: first, at T s ~ 
156 K, the finite local O domains begin to form side 
by side with major T domains in a slightly disordered 
manner; second, at T P ~ 148 K, the local O precursors 
have grown into long-range O domains. The tempera- 
ture range between T s and T P may be considered as a 
coexisting regime of the T and O phases. The split- 
ting of the (HKL)i (P4/nmin) reflection into twinned 
{H + K,H-K, L)/(H-K,H+K, L){Cmma) reflections 
is a sensitive measure of the T-to-0 structural transition. 
To track the O splitting, we performed a high-resolution 
synchrotron X-ray powder diffraction study with pulver- 
ized LaFeAsO single crystals. Part of the data are dis- 
played in Fig.Hwhich clearly shows the splitting of (220)t 
to (400)o/(040) and (322) T to (512) D /(152) D at 130 K 
below the O-T transition. It is noted that the (220) T 
reflection at 158 K is broader than the (400)o reflection 
at 130 K. This may indicate a distribution of d-spacings 
due to O-T fluctuations above T s . In Fig.[5j we show the 
temperature variation of the FWHM of the (004) / T and 
(220) T reflections obtained from a Lorentzian fit function 
to the Bragg peaks shown in Fig.|4j While the width vari- 
ation of the (004) / T is negligible within errors, the varia- 
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FIG. 6: (color online) (a) g-dependence of the collected fac- 
tor C B (q) (circles) in Eq. (2) at 9.5 K as described in the 
text. The nuclear Bragg reflections were chosen to cover the 
g-range of the measured AFM reflections as possible. The 
intensities of nuclear reflections were integrated from their 
respective rocking curves, taking into account the domain ef- 
fect, and normalized to the corresponding structure factors 
and the Lorentz factors. The collected factor Cs(q) contains 
the DWF and other effects such as extinctions and absorption. 
The solid line is the best DWF-like curve that best fits the 
data and is used to normalize the AFM peaks to obtain the 
magnetic form factor of iron in single-crystal LaFeAsO. (b) 
Measured magnetic form factor of iron (circles) at the AFM 
phase (9.5 K) of single-crystal LaFeAsO. The lines as indi- 
cated are the best fits of the data with the magnetic moment 
M (/xb) multiplied by the calculated form factors (jo) of Fe, 
Fe 2+ and Fe 3+ with spin contribution only |16| . and the form 
factor of a homogeneous sphere (details in text), respectively. 
The fitted results are shown in Tab. IIHI The deviation of A 
and B points from the smooth curves is discussed in the text. 



tion of (220)t is more significant showing a maximum at 
-150 K. The normalized FWHM of the (220) T reflection 
to that of the (004) D / T reflection displays an asymmetric 
temperature variation with a larger FWHM above the 
O-T transition [Fig. (5lb)]. This indicates remnant or- 
thorhombicity at temperatures as high as 300K. This is 



consistent with Ref. [T2] where a finite orthorhombicity 
remains visible above T P up to 200 K. Similar observa- 
tions have also been found in SrFe2As2 [HI [13]. We thus 
argue that both T and O phases may coexist dynami- 
cally due to the strong magnetic correlations in a certain 
temperature range [J2] above T N , which can be ascribed 
to a spin nematic phase [T7TU9| . We, therefore, associate 
T P = 148(1) K as the O-T phase transition temperature 
in the single-crystal LaFeAsO. We argue that this tran- 
sition temperature (T P ) may vary slightly, depending on 
the sample quality 19J, size and shape, and the cool- 
ing/warming protocols. 

Figures (3^b) and (d) show the temperature dependence 
of the rocking-curve integrated intensity of the AFM 
(103) reflection. The magnetic peak appears at ~140 K, 
indicating the formation of the AFM structure as shown 
inFig.Qb). 



B. The magnetic form factor of iron 

To obtain the magnetic form factor of iron in LaFeAsO 
we collected the integrated intensities of the rocking 
curves of nuclear and magnetic Bragg reflections to the 
highest order possible with the HB-lA spectrometer. The 
nuclear reflections, necessary to obtain the magnetic mo- 
ment size, were chosen to be able to cover as much of 
the g-range of the measured AFM reflections as possi- 
ble. This reduces any errors due to geometry corrections, 
Debye-Waller factor (DWF), absorption effects, and oth- 
ers. The integrated intensity of the rocking curve of a 
nuclear Bragg reflection at a reciprocal lattice vector q 
in a crystal is given by: 



4 = ^$oWI^(q)| : 

= C M (q)|F M (q)| 2 
sin(20) 



A 3 
2/xsin(20) 



,-2ff 



(1) 



where V is the scattering volume of the crystal, vo is the 
unit-cell volume, &o(0) is the beam flux at the angle 9, 
Fjs(q) is the structure factor, A is the beam wavelength, /j, 
is the absorption length, sin(20) is the Lorentz factor for 
a rotating crystal, and 2W = q 2 (uq) 2 is the DWF. We 
collect all the constants and unknown q-dependent fac- 
tors in Cn(<z), including the DWF. The integrated inten- 
sities of the rocking curves of the chosen nuclear Bragg 
reflections were normalized to the corresponding struc- 
ture factors and the Lorentz factors. The data were also 
corrected for the fact that the rocking curve of nuclear 
Bragg peaks in our measurements includes an equivalent 
twinned (/i00) and (QkO) domains [20], for normalization 
of the magnetic reflections which are due to the (/i00) 
domain only. The results are shown in Fig. J6^a) (circles) . 
We assume Ctt(q) is a DWF-like function (namely, a gaus- 
sian) and fit the data by the nonlinear square technique 



to obtain a smooth Cn(<z) function shown as a solid line 
in Fig. [^a). 

Similarly, the scattered intensity of the rocking curve 
of a magnetic Bragg reflection can be expressed as: 

h = C M (q)( 7n re^5)VM 2 (|q|)l^(/ifc0| 2 sin 2 a sI ^,(2) 

where C M (q) = C N (q)/4 due to the fact that the AFM 
unit cell is doubled along the crystallographic c axis in 
comparison with the nuclear one [Fig. [lib)], 7 n = -1.913, 
r e = 2.81794 x 10~ 5 A is the classical electron radius, 
/(i(|q|) is the magnetic form factor at the magnetic re- 
ciprocal lattice (q), \F K (hkl)\ = |£ e 2™(^+fc % +^)| = 8 
(j = 1 to 8) where (xj yj Zj) represents fractional co- 
ordinates of the jth atom in the AFM unit cell, and 
sin a = 1 — (q • /2) 2 where q and ft are the unit vec- 
tors along the scattering vector and the direction of the 
moment, respectively. Our attempt to calculate the mag- 
netic form factor with various directions of the ordered 
magnetic moment yielded an irregular and unphysical 
form factor, except when the moment points along the 
O a axis, namely, (1 — (100). This is strong confirma- 
tion that the AFM moment direction is along the O a 
axis. In Tab. |ll]we list all the observed magnetic reflec- 
tions, their integrated intensities, and the corresponding 
values of \F K (hkl)\ 2 sin 2 a. Using these parameters and 
Eq. (pi) we calculated the value of gSfn(\q\) as shown in 
Fig. pFb) (circles). To obtain the behavior of the form 
factor we fit the data in Fig. ffl(b) to different model form 
factors of iron in different environments. In general, we 
find that the form factor behavior is very similar to the 
one recently measured and calculated for SrFe2As2 |21j . 
Fitting the data with the theoretical form factors of Fe, 
Fe 2+ and Fe 3+ ((jo), spin contribution only) in Ref. [16], 
we obtain the corresponding average AFM moment sizes 
M (/Lt B ) as listed in Tab. |IIl| The experimental points in- 
dicated as A and B in Fig. 16Tb) deviate from the smooth 
curves, which may be due to systematic errors or due to 
some subtle features of the magnetic structure and its 
dynamics that are not captured by our structure factor 
dct crminat ion . 

As a simple, heuristic model, we assume that the 
magnetic moment is homogeneously distributed within 
a sphere of radius R, for which the form factor can be 



simply calculated by f(q) 



osin(g.R) 



-qR.cos(qR) 



!,,liy 



Fitting 



our measured data [circles in Fig. K3[b)] with this equation 
yields an average moment size of (0.74±0.04) fi B with a 
radius R = (0.63±0.04) A. This radius is almost the same 
as the effective ionic radius of Fc 2+ (0.63 A) 22 with a 
coordination number of 4, suggesting that the oxidation 
state of iron in LaFeAsO is closer to that of Fe 2+ . It is 
pointed out that the effective ionic radius depends on the 
particular electron configuration as well as the surround- 
ing ions (As) and the relative amount of ionic bonding. 
The average magnetic moment is ~0.8 /xb in the AFM 
phase of single-crystal LaFeAsO in disagreement with 



the reported smaller value of ~0.37 /U B (polycrystal) [BJ, 
but more or less near the recently reported ^0.63(1) 
Mb (polycrystal) [12]. In addition, the obtained mo- 
ment size in this study is almost the same as the first- 
principles-optimized 0.87 /ib in the magnetically frus- 
trated state [23], the ~0.83 /i B in an oxygen deficient 
LaFeAsO (polycrystal) [H], and the -0.8 fi B in CeFeAsO 
(polycrystal) [25] . First, it is pointed out that there is 
a large difference in the properties of polycrystalline and 
single-crystalline samples, especially for those containing 
easily-sublimating elements or oxides, e.g., the nominal 
Lao.875Sr .i25Mn03 [2"6"H2"5] . Second, single crystals are 
believed to be more stoichiometric [26] [28] and main- 
tain translational symmetry over macroscopic distances, 
in contrast with polycrystals, thereby providing more re- 
liable information on the structures and intrinsic prop- 
erties. The small magnetic moment size ~0.8 ^b (com- 
pared to Fe 2+ in an insulator, S = 2 with an average 
ordered localized moment 4 Mb) is comparable to that of 
the ^Fe 2 As 2 (A = Ca, Sr, Ba, "122") suggesting an ef- 
fective S « i on the iron sites [2T], and to some extent 
points to the itinerant character of the magnetism in this 
system. This puts the "1111" compound on a similar par 
with the other "122" parent compounds with a magnetic 
moment that is in the range of 0.8-1.1 /ib- This suggests a 
sharp jump in the magnetic moment of Fe ions as a func- 
tion of the Fe-As distance to about 1 ^ B , as predicted 
from first-principle calculations of bulk zinc-blende FeAs 
[29]. 

There are a few possible explanations to the small mag- 
netic moment measured by neutron diffraction techniques 
in iron-pnictides. The most accepted one is that it is due 
to the itinerant character of magnetism in this system. 
On the other hand, in a local moment picture this may 
be due to the dimensionality of the system or due to 
competing nearest neighbors (NN) and next-NN (NNN) 
exchange interactions. The Fe layers are weakly coupled 
making each layer behave as a quasi-2D system. The 
2D systems, especially those with isotropic NN coupling, 
exhibit strong magnetic fluctuations that can lower the 
quasi-static ordered magnetic moment. Second, the un- 
usual spin arrangement observed in iron-pnictides indi- 
cates the presence of strong competing and conflicting 
interactions between NN (Ji) and NNN (J 2 ) that can 
lead to strong magnetic fluctuations and a reduced static 
moment. 

To summarize, employing both neutron and syn- 
chrotron diffraction techniques to explore the details of 
the structural and magnetic properties of single-crystal 
LaFeAsO, we found: (1) The O-T structural transition 
occurs at T P ss 148 K, but the finite local O precursors 
appear to form already at T s s=s 156 K. We argue that the 
T and O phases may coexist in the temperature range of 
T s and T P , and at T P the long-range O phase has formed. 
(2) The AFM structure forms at T N w 140 K upon cool- 
ing with the iron moment direction along the crystallo- 



graphic a axis in the O structural regime. (3) The av- 
erage AFM moment size is comparable with that of the 
^Fe 2 As 2 (A = Ca, Sr, Ba), e.g., M = (0.82±0.03) [i B with 
a form factor of Fe 2+ with spin contribution only. This 
moment size is significantly larger than the previously 
reported values. More detailed studies of the form-factor 
will be required to determine the possibility of the re- 
cently predicted spin spatial anisotropy [21]. This study 
shows that to a large extent the properties of LaFeAsO 
are very similar to those of the so-called "122" system. 
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